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Abstract 



The elastic energy functional of a thin elastic rod or sheet is generalized 
to the case of an M-dimensional manifold in ^-dimensional space. We de- 
rive potentials for the stress field and curvatures and find the generalized von 
Karman equations for a manifold in elastic equilibrium. We perform a scal- 
ing analysis of an M — 1 dimensional ridge in an M = N — 1 dimensional 
manifold. A ridge of linear size X in a manifold with thickness h <C X has a 
width w ~ 

/ji/3_x-a/3 an d a total energy E ~ [ih M ' (X/h) M 5//3 , where \i is a 
stretching modulus. We also prove that the total bending energy of the ridge 
is exactly five times the total stretching energy. These results match those of 
A. Lobkovsky [Phys. Rev. E 53, 3750 (1996)] for the case of a bent plate in 
three dimensions. 
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I. INTRODUCTION 



The crumpling of a thin elastic sheet is mediated by the formation of a network of narrow 
ridges . Plastic deformation of the material in the neighborhood of these ridges leads to the 
ubiquitous linear scars in crushed paper, aluminum foil, and car bodies ^ was recently 

discovered by Witten, Lobkovsky, and others that this phenomenon can be accounted for 
using linear elasticity theory &-[7|. The scaling laws for a ridge were first derived by Witten 
and Li using a Flory type argument. A ridge of length X in a sheet of thickness h and 
Young's modulus Y was found to have a total elastic energy E ~ Yh 3 (X/h)^ 3 and a width 
w ~ /i 1//3 X 2 / 3 . Lobkovsky confirmed this result with a scaling analysis of the von Karman 
equations describing a thin, semi-infinite strip with a single ridge. He also verified these 
results with detailed simulations. 

One important result of this analysis is the discovery that the stresses and curvatures 
decay rapidly to zero in the direction transverse to the ridge. The length scale of this decay is 
the ridge width. We therefore interpret ridge formation as a confinement of the elastic stress 
field ||. Although there is a qualitative appreciation that confinement is the result of the 
competition between the in-plane strains and the curvatures of the plate, a deeper theoretical 
understanding is still lacking. In particular there is no proof of confinement under generic 
boundary conditions. We have been working towards this goal. In a companion paper we 
prove that a thin plate must have regions of nonzero strain if it is to fit into a small sphere 
|| . It remains to be shown that these regions necessarily assemble themselves into a network 
of ridges. 

To understand the causes and consequences of stress confinement, it is useful to examine 
the higher dimensional analogs of a crumpled sheet. For example, a thin plate crumpled in 
four dimensions doesn't have to stretch, and we expect there is no stress confinement in this 
case. In this paper we present the simplest field theory describing the strains and curvatures 
of a deformed elastic manifold in higher dimensions. We also perform a scaling analysis for 
a ridge in these systems. 
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There is considerable precedent for examining higher dimensional systems for insights 
into membrane elasticity. Most notable are studies of the so-called "crumpling transition" 
of a thin elastic membrane in thermal equilibrium fT0HT2|. That work focuses on the way 



thermal fluctuations and self-avoidance renormalize the elastic constants of a thin elastic 
sheet. The field theory is often developed in perturbation theory around a convenient, 
higher dimension. In this paper we work exclusively at zero temperature and the word 
"crumpling" refers to compression by external forces. 

In elasticity theory it is common to approximate a thin plate by its center surface, or 
centroid []13|,|l4[ . A three-dimensional plate is thereby described using a two-dimensional 
manifold. The elastic energy functional for the centroid is found by integrating out the 
components of the stress and strain tensors which are transverse to the long directions. The 
mathematical analysis of these approximations and their range of validity is the subject 
of the theory of thin elastic shells [I5]-|I7|. In this paper we use the standard methods of 



shell theory to derive the elastic energy of an M-dimensional manifold embedded in N- 
dimensional space. We treat the manifold as the centroid of an iV-dimensional elastic solid 
with an infinitesimal thickness h in N — M directions. The resulting energy functional 
has pieces quadratic in the strains, curvatures, and torsions of the manifold. We take a 
functional derivative of the energy to find the equations of static equilibrium. For a plate in 
three dimensions these equations are called the von Karman equations, first written down 



by Theodore von Karman in 1910 |1£|. We will refer to our general result by the same name. 
Our derivation has several new elements. 

It is noteworthy that the elastic energy of a thin plate may be written in terms of two 
scalar potentials. The stress function x, introduced by Airy in 1863, is the source of in-plane 
stresses []I9~| , |20"|| . The bending potential f is the source of curvatures. Similarly, studies of 
a deformed solid often use the tensor stress function Xa/3, introduced by Maxwell in 1870 



21| , |22[| . In this paper we present the generalization of these potentials for arbitrary M and 
N. Our derivation of the von Karman equations reveals a role for the stress function as the 
Lagrange multiplier of a geometric constraint in the energy functional. For M > 2 the stress 



function is a gauge field. 

Next we turn our attention to the ridge structure. The width of the ridge in a plate 
scales like w ~ 

^l/3 X 2/3 

so in the limit the ridge is approximately one-dimensional 

and straight. Analogously, the ridges in an M-dimensional manifold are expected to be 
approximately (M — l)-dimensional and to have no curvature. Indeed, simulations of a 
solid elastic ball crushed by a sphere in four dimensions show that the elastic energy is 
concentrated into flat, planar structures As mentioned above, Lobkovsky has done a 
thorough analysis of an isolated ridge for the case of a semi-infinite sheet in three dimensions 
0. We repeat his analysis for the case of an M-dimensional, semi-infinite manifold bent 
into a ridge in M + 1-dimensional space. We find that a ridge with linear size X has a width 
w ~ /i 1 / 3 X 2 / 3 and a total energy E ~ Yh N (X/h) M ~ 5 ^ 3 . Lastly, we prove that the total 
energy due to the curvature of the ridge is exactly five times the energy due to the strains. 
These results match correctly onto the solution for a bent plate in three dimensions |6,|23 



In Sec. II we review the differential geometry of a weakly strained M-dimensional man- 
ifold embedded in A-dimensional space. In Sec. Ill we derive the elastic energy of this 
manifold as the thin limit of an A-dimensional elastic solid. In Sec. IV we present the 
generalization of the stress function and the bending potential. Then we make a variational 
derivation of the von Karman equations. In Sec. V we generalize the scaling analysis of 
Lobkovsky to a ridge in M > 2. In Sec. VI we summarize our conclusions. 

II. DIFFERENTIAL GEOMETRY REVIEW 

In this section we review the differential geometry of a weakly strained M-dimensional 
manifold M. embedded in A-dimensional space 8l N . By weak strains we mean (1) the strains 
are small compared to unity and (2) the derivatives of the strains are small compared to the 
other relevant inverse lengths (curvatures and torsions). These are the usual assumptions 
of thin plate theory jTJ]. Note that they do not prohibit arbitrarily large deformations of 



the manifold. With these assumptions the math simplifies considerably. A treatment of the 
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topics in this section using the full apparatus of differential geometry may be found in Refs. 

US 



The manifold is flat in the absence of external forces, so it can be parameterized by the 
Euclidean coordinate patch {x = a^e, G M. for i G [1, M] ; Xi = for i G [M + 1,-AT]} 
where {e^} are the Euclidean basis vectors. We will refer to (x a ; a G [1, M]) as the manifold 
coordinate patch and denote it with Greek subscripts. Any deformation of M. can then be 
represented as a continuous map r(x a ) from the manifold coordinates to $l N . At each point 
on the deformed manifold there is an M-dimensional tangent space spanned by the tangent 
vectors t a = 6 a r. The metric on the manifold is then g a p = t a ■ tp and the strain tensor is 
u a p = (l/2)(g a/3 - 8 a p). 

We treat all relevant quantities to lowest order in the strains. This immediately gives us 
9ap — 8a/3+0(u), so there is no need to distinguish the covariant components of a tensor from 
the contravariant components. The Christoffel symbols are rig = dpu ai + d a u^ — d^Uap + 
0(u 2 ) P5 |. Thus, D a — d a + 0(u) and covariant derivatives are just partial derivatives to 
leading order. Geodesies are approximately straight lines in the manifold coordinates. 

The extrinsic curvature tensor for the manifold is defined K a p = D a tp pa d a tp. It is 
straightforward to show that the components of this tensor are normal to the tangent space. 
We start with 

t a ■ Kpy = t a ■ dptj 

= dp(t a ■ ty) — dfit a ■ try 

= 6/3(5 a y) — K a /3 ■ ty 

= —ty ■ K a /3 (1) 

This quantity is therefore odd under a cyclic permutation of the indices. Three such per- 
mutations gives t a ■ Kpy = —t a ■ Kpy = 0. 

We choose a set of orthonormal basis vectors {n^ a \x/3)} to span the (iV — M)-dimensional 
normal space at each point on the manifold. Note that we use Greek-in-parenthesis for 
the normal index (a) G [M + 1,7V]. In this basis the extrinsic curvature tensor becomes 



K a p = C^p where C^g = ■ K a /3 and summation over repeated indices is implied. We 
will refer to the N — M tensors C^I as the normal components of the extrinsic curvature 
tensor. Note that C^l = Cj^ and K a p = Kp a since K af} m d a dpr. 

It is useful to expand the derivatives of the normal vectors in the full basis {t a ,n^} 

da nW = _ C W r 7 - TfW ftM (2) 

where we have defined the torsions t^^(x) = —n^ ■ d a h^ and where we have used 
= d a {t-y ■ fi^') = t 7 • d a n^ + C^j . Eq. (|2|) is the generalization of the Weingarten map for 
a plate [p5f] . Taking one derivative of ■ fffl = 5( a )03) gives us the antisymmetry property 
for the torsions t^W = -t^. 

The last quantity we will need is the intrinsic curvature tensor R a pfi V — ■ Kp v — K av ■ 
Kfa p5] . In the normal basis, this becomes 

The intrinsic curvature tensor is related to the strain tensor via the generalization of Gauss' 
Theorema Egregium 

RaPuvM = -U afJ ,,l3u + U aU) p^ - Up u>ctfl + Uf3^ au + 0(u 2 ) (4) 

where indices to the right of a comma indicate partial derivatives. Eqs. (Q) and (|4]) together 
give R a pnv\P\ = Ra/3fiu[u}- This is one version of the geometric von Karman equation, so 
called because it expresses the geometric constraint relating the extrinsic curvature and the 
strain. It is straightforward to verify this equation by substituting the definition of the strain 
tensor into R a /3^y[u] and differentiating. 

The intrinsic curvature tensor will be most useful to us in the linear combination 

GaP = Rotvfiv — -^apRfivnv (5) 

This is the Einstein curvature tensor, familiar from general relativity [p6[] . It is symmetric 
G a p = Gf3 a and satisfies the conservation law d a G a p = 0. Taking the appropriate contrac- 
tions of the geometric von Karman equation gives us 
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G a p[C] — G a p[v\ 



(6) 



III. THE ELASTIC ENERGY FUNCTIONAL 

In this section we obtain an expression for the elastic energy of an M-dimensional man- 
ifold M. via the thin limit of an N- dimensional solid M . Versions of this calculation for a 



rod and plate may be found in Refs. |L3[ and [14 



To begin we consider the elastic energy functional for an arbitrary ^-dimensional solid. 
We keep the assumption of small strains used in the last section, but we relax the condi- 
tion on the derivatives. As before, there is a Euclidean coordinate patch (xj) covering the 
undeformed manifold {x = x^i G M for i G [1, N]}. This is the material coordinate patch, 
denoted by Latin indices. Under the application of external forces the solid assumes an 
embedding f(x). The tangent space of the solid is the full 3l N and the tangent vectors are 
U(x) = dir. The metric is g%j{x) = U ■ tj and the strain tensor is Uij = (1/2) — <%). 

Two consequences of u^- <C 1 are (1) the volume element of the deformed solid dr N (x) is 
well-approximated by the volume element of the undeformed solid dx N and (2) the elastic 
energy of the material only needs to be calculated to second order in the strains. The most 
general energy functional quadratic in the strains and consistent with an isotropic material 



is 



E[u] = /^(K- + ^4) (7) 



where \x and A are the Lame coefficients [|14[|. It is useful to rewrite this equation E 



(1/2) / dx N OijUij where 

(Tij(x) = 2/j.Uij + X6ijU kk (8) 

is the stress field conjugate to Uy. The stress field satisfies the conservation law D^a^ = 0. 

By analogy with the treatments of a rod and sheet, we assume that J\f — A4 x B N ~ M {h) 
where B N ~ M (h) is an (JV— M)-ball of infinitesimal radius h. The choice of a spherical "cross- 
section" is important to preserve the full rotational symmetry in the normal space of A4. 



The material coordinate patch becomes {x = x^i E M for i e [1,M] ; Yh=m+i x 1 ^ ^ 2 }- 
We refer to the long directions as the manifold coordinates and denote them with Greek 
indices. We refer to the short directions as the normal or transverse coordinates and denote 
them with Greek-in-parenthesis. For clarity we relabel the transverse coordinates = X( a ), 

SO that (Xi) = (Xi, X 2 , . . . X M , C(M+1), C(M+2), • • • C(JV))- 

The M- dimensional surface satisfying = is the centroid of A/". When the transverse 
degrees of freedom are integrated out it is the centroid which becomes the manifold M. 
Under the application of external forces the centroid deforms to some equilibrium embedding 
r°{x) = r(x, ( = 0). All the quantities discussed in the previous section are well-defined with 
respect to this embedding. With the exception of C^j and t^ x \ we will denote quantities 
calculated on the centroid with a superscript c. 

The first step in deriving the elastic energy functional for the centroid is to make a Taylor 
expansion of the embedding f(x, () in ( 

f(x, o = nx) + c w aP> (x) + IcooCm (*) + ••• (9) 

where a^\x) = d^f(x, C)lc=o an d b^ v \x) = d( p )d^ u )f(x, ()\(=o- Recall that in Sec. II we 
had the freedom to choose an arbitrary set of torsions due to the rotational symmetry of 
the normal space. Here we make the natural assignment #(1) = aP>/\aW\. With this 
identification the torsions of the normal basis are the torsions of the deformed solid. 

To calculate the energy we need to solve for and ftMM in terms of u^g, Cj^j , rjf^ u \ 
and their derivatives. We make the following assumptions: (1) u c af3 <C 1, Cj^j <C 1/h, and 
r MW <^ i//^ w e w iH see that these are necessary to satisfy the small strain condition 
Uij « 1. (2) The smallest length scale £ over which the strains, curvatures, and torsions 
vary satisfies I ^> h. We therefore write the most general expressions consistent with the 
rotational and reflection symmetries of the problem to first nontrivial order 

a (fl) = (l + a lU c aa )n^ 
gMM = 6l (Cfg n M + CM «W) + b 2 5^\C^ n (A) ) (10) 



8 



where ai, b\ } and 62 are dimensionless constants to be determined. The corrections are of 
relative order 0(u,hC,hT,andh/£). A more detailed account of the derivation of Eq. ( |i0|) 
may be found in Appendix A. Note that and feMM are independent of t a and t^^ v \ 
primarily because the torsions are antisymmetric under fi <-> v. If the cross-section of M 
isn't rotationally symmetric, then there are additional possibilities in Eq. (]T0|) . This may 
couple the torsions to the curvatures in a nontrivial way and complicate the resulting theory 
considerably. 

The tangent vectors to first nontrivial order are 

d a r = it + C W {-CS? r 7 - r a WW ftM} 
9 (a) r = (1 + aiu^) n w + &iCoo(Cg + C$ n< a >) + &2C(a)(C# n< A >) (11) 

and the components of the strain tensor are 

u am =-lCv>-ft m (12) 

From this it is clear that our assumption (1) above is equivalent to the small strain condition 
Uij <C 1. Also note that the transverse derivatives of the strain tensor d( a \Uij are of 0(C, r) 
and are not negligible. Only the manifold derivatives d a u^ can be safely neglected to leading 
order, as assumed in Sec. II. 

We derive the energy for a portion of the manifold far from the regions where external 
forces are applied. We therefore have the boundary condition cr( a )(p)\(;=h — 0. Combined 
with the conservation law = 0, we have ff[ a )(p) = everywhere. This condition specifies 
u (a)(f3) uniquely. Referring to Eq. (§) 

a (a)(l3) = 2/iU( a )(/3) + A5( a )Q3)1ijfejfc = (13) 

so M( Q )(/3) ~ <5(a)(/3) an d b\ + 6 2 = 0. Substituting the trace of the strain tensor 

u kk = u c aa - C M C<g + (N- M)( ai u c aa + hC M C^) (14) 
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into Eq. (|13D gives us a\ = —bi = — cq where Co = A/(2yU + (N — M)X). Thus 



-c 5 {am (u c 77 -C ilx) C^) (15) 
u kk = Cl (<4 7 - C(m)C^) (16) 

where c x = 2/V(2/x + (iV — M)A). 

Substituting Eqs. (|T2| ) and ( |16D into the expression for the strain Eq. (|j) gives 

The elastic energy Eq. (0) becomes 



+ c 4k 7 -C(,)^ ) ) 2 + ^(c ( ,)^^) 2 



(18) 



2 v 77 77 ' 1 ^ 

Since the transverse coordinates are being integrated over B N ~ M , terms odd in M ) vanish 
and we have finally 

£[«,C,t] = J M dx M {^(«^) 2 + c « a ) 2 ) 

where 



(19) 



V c = V / ^C 7V_M (20) 

JBK-M(h) 

is the effective stretching modulus and 

« = V I dC N - M Cl) (21) 

is the effective bending modulus of the thin manifold M.. The integral in Eq. (|20"D is just 
the volume of a sphere with radius h 

d( N-M = Kd^ (22) 

B d {h) a 
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where S d = 2n d ' 2 /T(d/2) is the area of a unit sphere in d dimensions. The integral in Eq. 



(0) is 



JB d (h) 



d/-1 
(1) 



d = l 
d = 2 



^h d + 2 B(3/2,d-2)S d _ 1 d>2 



where B(a, b) = T(a)T(b) /T(a + b) is the beta function [ pTfl . 
We can rewrite the elastic energy using conjugate fields 

E[u, C,t] = \I dx M {a^u aP + MgcM + 7^)M r WM } 



where 



77 1 



is the resultant strain field, 



M^(x) = 2K(C^ + c 6 aP C^) 



is the bending moment field, and 



(23) 



(24) 



(25) 



(26) 



(27) 

is the torsional moment field. Eq. (|24]) is the full elastic energy functional for a thin elastic 
manifold. We will frequently refer to the term quadratic in the strains as the stretching 
energy and the term quadratic in the curvatures as the bending energy. 

We henceforth drop the superscript c and assume that all quantities refer to the centroid 
manifold Ad. 



IV. THE POTENTIALS AND THE VON KARMAN EQUATIONS 



A. The case M > 2 



Note that Eq. does not explicitly couple the strains to the curvatures of the manifold. 
The strains and the curvatures are implicitly coupled because they are both defined via 
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derivatives of the embedding f(x). However, a naive functional derivative of Eq. (|24j) with 
respect to u a/ 3 gives the trivial and incorrect result a a p = 0. Previous authors, working with 
a thin plate in three dimensions, have solved this problem by working in a special coordinate 
system which is approximately tangent to one point on the manifold. In this frame, known 
as the Monge representation, the embedding is r(x) = [x a +u a (x),w(x)] and the derivatives 
of u(x) and w(x) are assumed to be small everywhere. To leading nontrivial order the strain 
tensor is u a p = (l/2)(u a ,p + up )(X + w^w^) and the extrinsic curvature tensor is C Q/ g = w }(X p. 
Functional derivatives are then taken with respect to u a and w and the correct equations 



are obtained 28 



We choose to work instead with the field variables u al 3 } C^) , and their potentials. The 
advantages are (1) we work exclusively in the manifold coordinates, so there is no need for 
an approximately tangent frame, (2) it is easier to treat the boundary conditions, and (3) 



we discover a new interpretation for the stress functions of Airy and Maxwell |19| , pT| , p2| . 
Due to some small differences, we focus here on the case M > 2 and return to the case 
M = 2 in the next subsection. The coupling between the strain and the curvature is 
completely accounted for by the geometric von Karman equation Eq. (|6|) . We therefore add 
the Lagrange multiplier term 

E x [u, C, X \ = \I dx M Xa p(G a p[C] - G a(3 [u}) (28) 

to the total elastic energy Eq. (|24]) . We will see that the Lagrange multiplier Xap{ x ) is the 
tensor stress function. 

One may ask why it is sufficient to use the Einstein curvature tensor G a p rather than 
the full intrinsic curvature tensor Rap^v G a p is symmetric in a (3 and constrained by the 
conservation law G a p j0l = 0, so a naive count of the independent degrees of freedom gives 
M(M — l)/2. R a pn U is symmetric in a ^ (3 and \i <-> v, and antisymmetric in (a/3) <-> (fJ,v), 
so in principle it has M(M — 1)(M 2 — M + 2)/8 independent components. However, Eq. 
(f|) shows that for small strains the intrinsic curvature tensor is linear in the strain tensor 
u a p. The strain tensor is symmetric and constrained by the conservation of the resultant 
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stress tensor, so it has MiM — 1)/2 independent components. The Einstein curvature tensor 
is therefore the most economical choice. The alternative forms for the Lagrange multiplier 
term Xa/J/if-Ra^ and XapRa^p^ both yield the correct von Karman equations for the fields 
a a p and C a X p , but the Lagrange multiplier is not identical to the stress function. 

If the normal basis has zero torsion we can define a bending potential for each normal 
component of the extrinsic curvature tensor. In the remainder of this paper we assume 
there are no external torsional moments acting on the manifold. Because the torsions are 



not coupled to the strains or to the curvatures in Eqs. (j24|) and (|2q) , the solution is simply 
= 0. With this we can prove the Codazzi-Mainardi relation 

= n< A > ■ {d a dpQ + (-(7$ t s ) ■ Kp, 



n (X) ■ (dpK, 



07 j 

= d p CSj (29) 

We have used the simplified Weingarten map d a n^ = —C^) tp in the second line and the 
orthogonality condition t a ■ h^> = in the third line. The Codazzi-Mainardi relation is 
analogous to the zero-curl condition on a vector field. It allows the definition of a scalar 
potential via C$ = d a dpf^ x l 

There is a novel form for the Einstein curvature tensor which greatly simplifies the 
variational derivatives taken below. We begin by defining the double curl operator, valid for 
M > 2, 

[d-C-JaPnv — Jj^j- _ ~y € a -y flTl ... TM _ 3 €pSi'T 1 —T M _ 3 d'yds 

+S tw dpd ll + Sp^dadv - S^dadp (30) 
where e T1 ... TM is the Levi-Civita tensor 
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+ 1 (n, r 2 , . . . , TAf) even permutation of (1, 2, 3, ... , M) 
— 1 (n, r 2 , . . . , tm ) odd permutation of (1, 2, ... , M) (31) 
otherwise 

The double curl is antisymmetric in a «-> /i and /3 v and symmetric in (a/x) (/3z/) 
(compare to R a /3^ u ). It satisfies d a {d.c.) a p^y = by construction. The Einstein curvature 
tensor may be written G[u} aP = {d.c.)^^^ or G[f] a p = (d.c.) a/3flv [(l/2)f^f^]. These 
expressions are easily verified by substitution. 

Now we can write the full expression for the energy functional, including the Lagrange 
multiplier 

E[u, /, X ] = J m dx Ml - {a a p[u}u a p + «(/$/$ + 

+ X Q/3 (rf.c.) Q ^(i/W/W - UfiV ) (32) 

Taking a functional variation with respect to u a p and integrating the Lagrange multiplier 
term twice by parts gives 



(33) 



5E = dx M bu^ip^ - {d.c.)^ ua pXap) 
Jm 

where we have used the symmetry of u a /3 and Xap- The first term in Eq. (|3^) gives a 
conservation law for the resultant stress tensor. Taking 5E/5u a p = 0, 

(J aj 3 = {d.C^p^Xtw (34) 

This is a restatement of the conservation law d a a a p = 0. Eq. (|34|) is the defining equation 



for Maxwell's stress function in M = 3 pi, 22]. We see that the stress function is a Lagrange 



multiplier. This interpretation persists even when M = N = 3 and the extrinsic curvature 
tensor is identically zero [take = in Eq. (ff^)]. Eq. ( |34"D also provides a natural 
generalization of the stress function to higher dimensions. 

One can verify by substitution that the stress tensor remains unchanged under the local 
gauge transformations Xa/3 — > Xa/3 + (l/2)(£ aj( g + £p, a ) where £ a (x) is an arbitrary vector 

14 



field [29|. The tensor stress function is therefore a gauge field with M(M — l)/2 physical 
degrees of freedom. This agrees with the fact that the stress tensor itself has M(M — l)/2 
independent components. 

The second term in Eq. ([33]) is a perfect differential. Using Gauss' Law to rewrite it as 
an integral over the M — 1 dimensional boundary of the manifold DM. gives 

5E\qm = — ■^y e aJf j iT1 ... TM _ 3 €pSuTi-TM-3 J dM 1 177>~{(8u^ u Xaf3,8 — SUfivjXap) (35) 

where rh is the unit outward normal defined in the tangent space of Ai. The application of 
this term to a specific problem depends on the boundary conditions imposed at dAi. 

Taking the functional variation of E with respect to and integrating by parts gives 

5E = [ dx M 6fW ( 2 «(l + c ) V 4 / (A) - f% (d-c)^xj 

J M 

+2n (SfWfW + Co W/M - (1 + Co ) tf/W/W,) (36) 

(M-3)! e a07-tri— TM-z e p5vrx — tm-3 ^7 Xa/3 ) 

where there is no sum on (A). The second and third terms are perfect differentials and may 
be written as a condition on dAi. Taking 5E/8f^' = 0, the first term gives 

2k(1 + co) V 4 / (A) = /S (d-c W X/- (37) 
which may be rewritten in the more familiar form 

M$U = (38) 

This is the force von-Karman equation, which expresses the balance of forces on a transverse 
section of the thin manifold. 

To complete our discussion of the von Karman equations, we write the Einstein curvature 
tensor in terms of the stress function. First invert Eq. ( pSj ) to get 

where c 2 = A/ (2fi + NX). Then 
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Ga/jfx] — — {(d.C.jap^ld.C.) fiuarXar ~ c 2(d-C.) a p^(d.C.) waT XaT} 
= ^ {(C3 - 1) (SapV^Xaa ~ Sa^Xar^T ~ ^X^afs) 

where c 3 = (M — 2)[2/i + (iV - M + 2)A]/(2/x + ATA). 

We thus have several alternative expressions for the von Karman equations, depending 
on which fields are most convenient. In terms of the extrinsic curvature and strain tensors, 
we have Eqs. (§) and (|38|) . In terms of the bending potentials and the stress function we 
have G aP [f] = G afS [x] and Eq. @7|). 



B. The case M = 2 

In this section we rederive the von Karman equations for a thin plate. Although the 
equations of the three-dimensional problem have been discussed in detail by several authors, 
our variational derivation of the force von Karman equation is particularly transparent 
jjl , 14 , 2(| 6[] . The only change from the previous section is that the tensor double curl 



operator is not defined for M = 2. Instead we use the scalar operator 

{d.c.) aj3 = e^e^d^dy = 5 a/3 V 2 - d a d^ (41) 

When M = 2, the intrinsic curvature tensor has only one independent component, which 
we take as the generalization of the Gaussian curvature 

kg[C] = -jR a pap[C} = C^C^ — C^C^ (42) 

It is straightforward to verify that Kg[u] = — {d.c.) a pu a p and Ka[f] = — {d.c.) a p (hf^fffl^)- 
The elastic energy is 

E[u, C, X ] = I dx M \ {<T a p[u}u a p + M$C$} + x(-kg[C) + k g [u)) (43) 

E[u, f, x] = J m dx M \ {a a ?iu}u a p + n{f%f% + cof%f%)} 

'1 
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M 2 , • ^ (44) 

+x{d.cU(-f$ff-u aP 



where the Lagrange multiplier x( x ) is a scalar field. Taking the functional derivative 
SE/5u a/ 3 = gives <r a/ g = {d.c.) a pX- Thus d a a a ^ = and x is the scalar stress function of 
Airy 

In terms of the stress function kq[x] — [(1 — c 2)/2/i]V 4 x and the geometric von Karman 
equation is 

hl-c 2 )V\ = -{d.c.U{f^ff) (45) 
The force von Karman equation is found via SE/Sfw = to be 

2*(1 + c ) V 4 / (A) = f% (d.c.U X (46) 
which reproduces Eq. ( |38f) . 

V. THE RIDGE 
A. Boundary Conditions 

In this section we discuss the picture of a ridge as a boundary layer and find a simple 
boundary condition which yields a ridge for general M. Previous analytic studies of stress 
confinement have been limited to the case of a thin plate (M = 2) in 3? 3 . In Ref. ||, 
Lobkovsky treated the case of a single ridge in isolation and analyzed the resulting von 
Karman equations to lowest order in the thickness. He found that a ridge of length X has 
an elastic energy E ~ / u/i 2 (X//i) 1 / 3 and a width w ~ /i 1 / 3 A 2 / 3 '. 

We develop the concept of the ridge as a boundary layer by treating the thickness h as a 
tunable parameter. We start with h = 0. Since the width of a ridge scales like h 1 ' 3 , the zero 
thickness limit of a ridge is a straight line of zero width. The geometry of the ridge, shown 
in Fig. 1, is two flat plates which meet at a nonzero angle D. The curvatures are obviously 
singular on this line and zero elsewhere. The intuitive reason for this behavior is that 
k/(i ~ h 2 , so there is no energy cost for curvatures when h = 0. When h is made nonzero, 
the plate can achieve a lower total energy by smoothing out the singularity and trading 
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stretching energy for bending. The resulting balance generates a new length scale, which is 
the width of the ridge. This picture will remain essentially unchanged when M > 2. The 
important point is that the ridge is a boundary layer which regularizes the h = singularity. 

Lobkovsky began with the semi-infinite strip {— oo < X\ < oo, —X/2 < X2 < X/2}. Then 
he assumed the presence of (unspecified) normal forces acting at the boundary x 2 = ±X/2 
sufficient to deform the strip into the ridge shown in Fig. 1. The boundary conditions are 
C a p(xi, =t^/2) = 0, a a p(xx, ±X/2) = 0, and /(xi, ±X/2) = a\x\\. The first two conditions 
are chosen for convenience. It is the third condition which determines the shape of the ridge. 
The potential / plays the role of the normal coordinate. The dihedral angle of the resulting 
ridge is D = n — 2a. 

The generalization of this geometry to M > 2 is straightforward. We limit our discussion 
to the hypersurface N = M+l, since this captures the most important features of the general 
case [J]. We take for our h = ridge the singular boundary between two M-dimensional 
regions with zero curvature. As discussed in Ref. |J, such a boundary can have no curvature 
in the material coordinates or in 3? . Our boundary conditions must be consistent with this 
ridge. 

We take for our undeformed manifold the semi-infinite domain {— oo < x\ < oo , x g e„ G 
C for a e [2, M]} (see Fig. 2). The material coordinate x\ is perpendicular to the ridge and 
C is an arbitrary, simply connected cross-section. The ridge is imposed by some (unspecified) 
normal forces sufficient to create the "kinked" potential f\dM = ot\x\\. For simplicity we 
take C a p\QM = and cr a/3 \g M = 0. The solution to the von Karman equations when h = 
is the singular ridge f(x) = a\x\\. 

B. Scaling analysis 

To begin a scaling analysis of the von Karman equations we need to decide which fields 
to work with. Since the stress function Xa/3 is a gauge field in M > 2, it is convenient to 
work directly with the stress tensor cr Q/ g. We use the bending potential p 1 ' instead of the 



18 



curvature tensor C^g because of the obvious advantages of a scalar. Since N = M + 1 we 
can drop the normal index. The geometric and force von Karman equations are then 

(d.c.) a/3lxu (f^f 7U ) = -{d.c.) a ^ u (a^ v - c 2 ^cr 77 ) (47) 

2«(1 + C )V 4 / = <Tapf,a(} (48) 

Now convert to the dimensionless quantities 

X 2 

f = f/X , a af3 = -a al3 and x a = x a /X (49) 

Zk{1 + c ) 

where X is a length scale characterizing the cross-section C. Including the conservation law, 
the dimensionless equations are 

d a a af3 = (50) 
{d.c.) af3flu (f^f\ u ) = e 2 (d.c.) af3flu (a^ - c 2 <^,<7 77 ) (51) 
V 4 / = a afi f >al3 (52) 

where e 2 = 2k(1 + co)/(/uX 2 ) ~ (h/X) 2 . As for a bent plate in three dimensions, we expect 
the ridge solution to be valid when < e < 1. Note that a naive count of Eqs. (|50|)-(|52|) 
gives M(M + l)/2 + N constraints on M(M + l)/2 + N — M field variables. This set is 
not overdetermined because the argument of the double curl has the local gauge freedom 
Onf — > Qnv + (1/2) (Cfi,v + £v,n) where £(x) is an arbitrary vector field. The geometric von 
Karman equation therefore represents only M(M — l)/2 independent constraints. 

The e = limit is the h = limit. As discussed in the previous section, the solution to 
this reduced problem is the singular ridge / = a\xi\. We might hope to solve for the e > 
ridge as a perturbation about this solution. However e multiplies the highest derivative of 
a a p. It is therefore a singular perturbation and naive approaches fail. 

To find the exponents which characterize the ridge, we make the following rescalings 

a n = e d a n , a m = e s cr m , = e*a^ 

/ = e^f , X\ = e^Xi and x a = e°Xa (53) 
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where we distinguish the coordinates parallel to the ridge with a barred Greek index a G 
[2, M}. Note that / ~ X\ is required by the boundary condition, so they scale with the same 
power of e. 

To solve for the exponents we start with the conservation law Eq. fl50|). Grouping terms 
of like order in e gives 

e^dxan + e^d^i = (54) 
eP-'hd^ + e-*9 a cr^ = (55) 

Assuming none of the rescaled quantities vanish, this implies s = 5 — (3 and t = 5 — 2(3. 
The rescaled force von Karman equation is 

(e 2 ^ + e°dirf = e^\a a ,l aP ) (56) 

If (3 > this equation is dominated by the lowest order terms in e@ as e — > and = 2(3 — 5. 
If (3 < then 4(3 = 2(3- 5. 

The rescaled geometric von Karman equation has M(M + l)/2 components, one for each 
component of G a p. For the scaling analysis they can be grouped into four classes Gn, Gm, 
Gaa (no sum), and G^s. Assuming none of the relevant terms vanish, we only need to 
consider one example from each class. The Gn component is 

e- 2P {l & pl & p - laafffl) = -e 2 - S {e 2/3 a^ + e 2P Ci a mM + e°c 2 (M - 2)a nM } (57) 

where C4 = 1 — (M — 2)c2- We have made extensive use of the conservation law a a p )a = and 
the symmetries of the problem to simplify the expression. The other components are derived 
in Appendix B. A careful analysis shows that all components yield the same constraint 
on the exponents. If this were not the case, the equations for the exponents would be 
overdetermined. If (3 > then Eq. (0) gives -2(3 = 2-5. If (3 < then -2(3 = 2-5 + 2(3. 
The constraints on the exponents are only solvable if (3 < 0. The unique solution is 

P = ~l> 5 = +1 ( 58 ) 
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These are identical to the exponents found by Lobkovsky for the ridge in (M = 2, N = 3) 
||. There are a variety of geometric intuitions associated with these exponents. The most 
important is that the ridge width is characterized by the transverse curvature via w ~ 1/ C\\ . 
Thus 



w 



e~ X ~ h^X 2 ' 3 (59) 



We refer the reader to Ref. || for a fuller discussion of the ridge geometry. 
It is instructive to consider the rescaled ridge energy to leading order in e 



E = / (X M e~?dx } 
Jm . 



1Z.M 
\yy t u 

I SI 

1 [2k(1 + C ) : , 2S \ ,_ 2 ~2 \ , .. ( e 



2 / «\2 ) (60) 



Gathering terms and using the geometric von Karman constraint 5 — 4/? = 2 gives 

£ = fih M (c s e- M -^ + c b e- M+2+ ?) (61) 

where c s and c& are dimensionless constants due to the stretching and bending energy re- 
spectively. In this form it is clear how the value of /3, and hence the width, is generated via 
the balance between strains and curvatures. As the width w ~ is increased, the bending 
energy decreases and the stretching energy increases. The ridge chooses the value (3 = —1/3 
which minimizes the total energy. Thus 

E ~ fih M e- M+5 ^ 3 ~ iih M {X/h) M -^ (62) 

Furthermore, Eq. fl6l"l) fixes exactly the ratio Cb/c s . At the minimum, 

= /i/i M ln( e ) (-5c s e- M ~ 5(3 + c b e- M+2+f3 ) _ =0 (63) 



BE 



(3=-l/3 



/3=-l/3 

This reduces to c&/c s = 5, which means the bending energy is exactly five times the stretching 
energy in an asymptotic ridge. As noted in the introduction, this ratio also holds for (M = 
2,N = 3) p3|. 
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VI. CONCLUSIONS 



In this paper we derive the equations of static equilibrium for an M-dimensional elastic 
manifold embedded in iV-dimensional space. We define the potentials Xap and on the 
manifold. These are the higher- dimensional analogs of the stress function x and the bending 
potential / of a thin plate in three dimensions. We find a novel interpretation for the stress 
function as the Lagrange multiplier of the geometric von Karman equation in the elastic 
energy functional of the manifold. 

We go on to consider the properties of an M — 1 dimensional ridge in an M > 2 di- 
mensional manifold. The scaling is essentially identical to that found by Lobkovsky for a 
ridge in M = 2 || . We find that a ridge of linear size X in a manifold of thickness h has a 
width w ~ /i 1 / 3 X 2//3 and a total elastic energy E ~ fih M (X/h) AI ~ 5 ^ 3 , where fj, is a stretching 
modulus. The scaling analysis also fixes exactly the ratio of bending energy to the stretch- 
ing energy in a ridge E bend / E stretch = 5. These results are valid in the thin limit /iCl 
Although our calculations are explicitly for a bent hypersurface N — M + 1, unpublished 
theory and simulations lead us to expect no change in the ridge exponents when N > M + l. 

The purpose of this work was primarily as an aid to future studies of crumpling in high 
dimensional systems. In particular, in future papers we will discuss the phenomenon of 
spontaneous ridge formation as a mechanism of stress confinement The elastic energy 
and ridge exponents derived here are an essential foundation for that work. We should point 
out that the scaling analysis in this paper is in no way a guarantee that these ridges will form 
in a crumpled manifold. The question of how the elastic energy is distributed is best resolved 
in combination with computer simulations and, for (M = 2, N = 3), experiments. As an 
example of the way ridge formation can fail, we note that when N > 2M a manifold with 
a free boundary can make its stretching energy zero everywhere. Since the ridge structure 
depends on the competition between bending and stretching energy, no ridge formation is 
possible. 
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APPENDIX A: COMMENTS ON THE TAYLOR EXPANSION 

Begin by expanding and in the full basis {t a: n( a )} 

gfcOM = jMM 4 + ft WM(A) ^(A) (A1) 

The coefficients ( a M aWM^WW^WMW) are f unc tions of u c af3 , eg, r^M, and their 
derivatives. We require that the expressions for the coefficients in terms of these quantities 
have the correct number of free manifold and normal indices (there may be an arbitrary 
number of contracted indices). This is necessary and sufficient for their correct behavior 
under reflections and rotations of the manifold and normal coordinates. 

Note that ft^M and ftMMW must be even under \i <-> v since b^^ v \x) = 
d(jj.)d(v)r(x, C)lc=o- This is why b 3 r^^ t a is not a valid term. Note also that = 0. 

We require that the coefficients have the correct units. Consider the following additions 
to the expression for 

a'^ = (o^gcg + a 3 rif)WrWM) ftM + for™^ + a 5 C^) f„ (A2) 

Since a*"> is dimensionless the constants dj must have units of (length) 2 . The only length 
scale available is the thickness h, so we write dj = djh 2 . The expressions are therefore 
of 0(h 2 C 2 , h 2 r 2 , h 2 Cr, and h 2 C/£) respectively. Although these terms are assumed small, 
they are not necessarily negligible compared to u c al3 . One can verify, however, that their 
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contribution to the energy is negligible compared to the curvature and torsion terms in Eq. 
Similar arguments lead to the form for ftMC"). 
In the theory of thin shells it is known that the centroid deformations alone are not 
sufficient to describe the behavior near the boundary of the shell. The full three-dimensional 
problem must be solved there. As a consequence, any energy functional derived via a Taylor 
series expansion in the thickness is not uniformly convergent near the boundary. For a 
detailed discussion of these considerations, we refer the reader to Ref. 0] and the references 
therein. 

APPENDIX B: SCALING OF THE FORCE VON KARMAN EQUATION 

Assuming none of the relevant terms vanish, we only need to consider one example from 
each of the three classes G m , G^a (no sum), and G & s. Dropping the tilde notation, the Gi 2 
equation is 



(Bl) 



2e P{f,V2.f,aa — f,mf,2a) 

= -e 2 ~\e w 0-12,11 + e P (T 12 ,aa + e /3 c 4 o-n,i2 + e^c^a^n) 
Defining dotted Greek indices a G [3, M], the G 22 equation is 

= -6^{c 4 (6 2 ^ lllll + e°o uM + + e^a^) (B2) 

+e 4/3 0-22,11 + e 2/ V 22 ,aa} 

and the G 2 3 equation is 

2e°(/,23/,ll - /,2l/,3l) + 2e~ 2/3 (/,23,/, / 3/3 ~ f,2$f,3$) moX 

(B6) 

= -e 2 - 5 (e c 4 (Tii,23 + e 2/3 c 4 a-^,23 + e 4l3 a 23 ,n + e 2l3 a 23 ,aa) 

If /3 > 0, then all three equations give the constraint —2/3 = 2 — 5 in agreement with the 
Gn component. If j3 < 0, then all three give = 2 — 5 + 4/3, again in agreement with G u . 
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FIGURES 




FIG. 1. The ridge in a semi-infinite strip when h = and h > 0. We have labeled the h = 
ridge to show the ridge length X, the dihedral angle D, and the manifold coordinate system (x±, X2). 




FIG. 2. The manifold coordinate system (xi,X2,X3) for the semi-infinite rod 3ft 1 x C. The 
location of the planar ridge is indicated in grey. 
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